应用Delaunay图映射与FFD技术的层流翼型气动优化设计

采用非均匀有理B样条(NURBS)基函数属性建立了任意空间的自由式变形(FFD)翼型参数化方法,进一步结合基于Delaunay图映射技术建立了结构对接网格变形模式,通过粒子群优化(PSO)算法进行参数化方法、网格变形模式以及计算流体力学(CFD)数值模拟技术之间的整合,研究、构建了气动优化设计系统,并对某型层流理念设计的高空长航时(HALE)飞机基本翼型进行气动优化设计。气动特性目标函数评估方法中,边界层转捩数值模拟技术采用γ-Re_θt转捩模型耦合剪切应力输运(SST)模式湍流模型。优化设计后翼型气动特性表明:采用相关技术建立的层流翼型气动优化设计系统对于层流理念设计的HALE飞机翼型的设计具备较高的优化效率。     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

某涡轮级三维流场非定常数值研究

基于耦合Lantry-Menter转捩模型的SST(Shear stress transport)湍流模型,对某涡轮级在3个雷诺数工况下的全三维流场分别进行了定常、非定常数值模拟,目的在于对比分析两种计算方法所得结果之间的差异。计算结果表明,进口雷诺数为4×104时和5×105时,定常与非定常计算的涡轮级效率的差异较大,最大相差0.6个百分点,而当雷诺数为8×104时,两者差异很小;非定常计算能够较好的模拟叶尖泄漏涡的发展过程、低雷诺数下上游尾迹与下游转子叶片吸力面分离边界层干涉诱导卷起涡的形成、输运及其所导致的叶片吸力面表面压力大幅波动等非定常流动现象。最后,给出了转子出口截面上周向与径向的非定常压力脉动的分布,为基于Lighthill声类比方法的计算气动声学研究奠定了基础。     

二维翼型微吸吹气减阻控制新技术数值研究

结合抽吸气转捩控制和微吹气湍流减阻控制的特点,探索了一种新的吸吹气减阻控制技术。使用Fluent求解器,并利用用户自定义函数(UDF)二次开发对其自带的Wilcox转捩模式进行了修正。在此基础上,数值研究了吸吹气控制对翼型阻力性能的影响。结果表明：在一定的吸气量范围内,吸气、吸吹气控制都能使翼型总阻力减小,且在同一雷诺数下,吸气控制能使翼型总阻力减小约3%,而吸吹气联合控制使翼型总阻力减小约16%。由此可见,吸吹气控制技术是一种行之有效的减阻控制技术。     

类X-51A飞行器非定常湍流精细模拟

针对类X-51A飞行器在超声速大迎角状态下存在的大范围非定常分离流动,开展了精细化湍流数值模拟研究。计算基于高阶格式下的延迟分离涡模拟方法（DDES）,来流马赫数为2.5,迎角为10°。分析了该复杂流场中存在的分离流动现象、分离流动诱导的气动特性变化规律以及压力脉动特点;其中重点研究了壁面压力脉动强度分布情况和监测点压力脉动频谱特性。分析结果表明:飞行器大迎角飞行时从侧缘诱导出明显的分离涡,并对尾部舵面产生干扰;受干扰尾舵表现出明显的非线性及非定常气动特性;分离涡的存在导致飞行器尾舵前缘等位置的壁面压力脉动显著增强,200~300 Hz的低频高幅值脉动可能会导致结构破坏。      

航空发动机压气机内气液非平衡冷却特性研究

针对航空涡轮发动机来流雾化冷却对压气机内气动脉动的影响,考虑壁面液膜成形和运动,基于欧拉-拉格朗日多相流方法解析气液两相热质非平衡传输过程,应用快速傅里叶变换方法将压气机性能参数随旋转周期演变规律的时域脉动敏感性转化为频域功率谱密度的直观分析。结果表明,压气机内气液非平衡蒸发相变易诱发流场在时间和空间上非定常的气动脉动,雾化冷却参数与总温比成线性关系,而与总压比和效率均成非线性关系。在雾化量0.5%~5%和雾化平均粒径1~9 μm范围内,较低的雾化量或较小的雾化平均粒径时,时域总压比的脉动程度更大;在较低的雾化量或较大雾化平均粒径时,时域总温比的脉动程度更强;而在较高的雾化量或较大的雾化平均粒径时,时域效率脉动程度更高。同时,雾化冷却量变化对湿压缩过程中流场的时域脉动敏感性程度大于雾化粒径变化。     

甲烷摩尔分数及初始压力对甲烷/RP-3航空煤油混合燃料燃烧特性的影响

采用定容燃烧实验装置获得初始温度为450K、初始压力为0.1~0.3MPa、当量比为0.7~1.5以及甲烷摩尔分数为0~0.8工况下甲烷/RP-3航空煤油混合燃料火焰发展特性图片、马克斯坦长度和层流燃烧速度等燃烧特性,分析甲烷摩尔分数及初始压力对甲烷/RP-3航空煤油混合燃料燃烧稳定性及层流燃烧速度的影响。结果表明:当量比为1.3时,随着甲烷摩尔分数增加,甲烷/RP-3航空煤油混合燃料燃烧趋于稳定,初始压力对燃烧稳定性影响较大,随着初始压力增加,燃烧稳定性变差。混合燃料马克斯坦长度随当量比增加而减小,当甲烷摩尔分数增加时,混合燃料马克斯坦长度减小趋势变缓,当初始压力增加时,混合燃料马克斯坦长度减小趋势明显变缓。混合燃料层流燃烧速度随当量比增加呈现先增大后减小的变化趋势。当甲烷摩尔分数为0、0.4和0.6时,随着甲烷摩尔分数增加,混合燃料层流燃烧速度逐渐增大,当初始压力为0.1、0.2、0.3MPa时,随着初始压力增加,混合燃料层流燃烧速度显著降低,随着甲烷摩尔分数和初始压力的增加,混合燃料层流燃烧速度峰值有向当量比大的区移动的趋势。      

测压孔径对下沉安装压力传感器测量的影响

对于高超声速飞行试验,飞行器在飞行过程中产生的气动加热会使其表面安装的传感器烧坏并导致试验测量失败。传感器下沉安装可以避免其与飞行器周围的高温气体直接接触,减轻传感器的热载荷。采用下沉安装传感器的方式对高超声速边界层的脉动压力进行了试验测量,研究了测压孔径对脉动压力特性的影响。结果表明,随着孔径增大,脉动压力强度减小,流场相关性增强。传感器下沉安装会引起空腔流动,随着孔径的增大,空腔流动对脉动压力测量的影响降低。在频带10~20 kHz范围内,测压孔径仅仅影响脉动压力强度,在频带1~10 kHz范围内,测压孔径对脉动压力波形周期的影响更为显著。非线性耦合将能量由高频小尺度向低频大尺度传递,导致发生自相互作用和非线性相位耦合的流动结构趋向低频,随着孔径增大,非线性相位耦合消失。此外,孔径越大,不同尺度结构的能量分布与齐平安装的相似性越高。     

不同导叶数下液力透平蜗壳内压力脉动计算

当离心泵作液力透平运行时其内存在振动现象,为了使液力透平能够稳定运行,在液力透平蜗壳出口设计不同数量的导叶,然后在不同导叶数下利用Pro/e软件建立几何模型,并利用ANSYS-CFD软件通过在蜗壳内沿周向和径向设置监测点,计算在不同导叶数下液力透平蜗壳内的压力脉动幅值,之后通过快速傅里叶变换将压力脉动的计算结果进行数据处理,分析不同导叶数下液力透平蜗壳内的压力脉动频域分布和液力透平机组内的振动情况.研究结果表明:当导叶数等于9时,蜗壳内周向和径向各监测点处的压力脉动主频幅值最小;蜗壳内径向各监测点处的压力脉动主频幅值和最大脉动幅值随着流量的增加而增大,但随着导叶数的增加其增大程度逐渐减小;当液力透平蜗壳出口添加导叶数为9的导叶时有效降低了液力透平机组内的振动和噪声,提高了液力透平机组运行时的稳定性.      

喷注壁面楔状流层流边界层速度与温度相似解

建立了壁面有喷注的楔状流层流边界层冷却数学模型,求解经相似变换得到的描述无量纲流函数和无量纲温度的常微分方程,获得了楔状流层流边界层无量纲速度和温度的相似解,给出了考虑壁面喷注边界条件的形式简洁、物理意义明确的无量纲流函数拟合式;用Runge-Kutta法求解无量纲速度与温度的常微分方程,获得了壁面有喷注的楔状流层流边界层无量纲速度和温度随楔形角、吹风参数、冷却介质温度的变化规律.通过计算0°,18°和36°楔形角的楔状流层流边界冷却速度与温度分布得到以下结论:楔状流速度边界层和温度边界层厚度都随楔形角的增大而变小,随着吹风参数的增大而增大;当冷却介质温度越低、吹风参数越小、楔形角越大时,靠近壁面处楔状流层流边界层内温度梯度越大.